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We propose two models of the evolution of a pair of compet- 
ing populations. Both are lattice based. The first is a compromise 
between fully spatial models, which do not appear amenable to an- 
' alytic results, and interacting particle system models, which do not, 

, at present, incorporate all of the competitive strategies that a pop- 

ulation might adopt. The second is a simplification of the first, in 
which competition is only supposed to act within lattice sites and 
■ the total population size within each lattice point is a constant. In 

a special case, this second model is dual to a branching annihilating 
random walk. For each model, using a comparison with oriented per- 
colation, we show that for certain parameter values, both populations 
will coexist for all time with positive probability. As a corollary, we 
deduce survival for all time of branching annihilating random walk 
for sufficiently large branching rates. We also present a number of 
conjectures relating to the role of space in the survival probabilities 
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1. Introduction. Natural populations interact with one another and with 
their environment in complex ways. No mathematical model can possibly in- 
corporate all such interactions and still remain analytically tractable. As a 
result, in order to understand the effects of a feature of a population's dy- 
namics, it is often useful to study "toy models." In this paper, we investigate 
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two such toy models that aim to parody the evolution of two populations 
that are distributed in space and competing for the same resource. Both 
of our models can be viewed as a compromise between fully spatial models 
which do not appear to be amenable to a rigorous mathematical analysis and 
interacting particle system models which do not, at present, incorporate all 
of the competitive strategies that a population of, say, plants might adopt. 

Although lattice based, our first model is highly reminiscent of the models 
in continuous space studied by Bolker and Pacala [2] and Murrell and Law 
[12], while admitting a rigorous mathematical analysis. It comprises a system 
of interacting diffusions, indexed by Z*^, driven by independent Feller noises 
and coupled through a drift term that reflects migration and competition 
(both within and between species). Our second model is much simpler: first, 
we suppose that the parameters governing migration of individuals within 
the two populations are the same and that competition between the popu- 
lations acts only within individual lattice sites; second, we suppose that the 
total population size within each lattice site is a fixed constant. When we 
further restrict to the symmetric case, in which the parameters governing 
the evolution of the two populations are the same, we shall exhibit a duality 
between this second process and a branching annihilating random walk. The 
latter is a process that has received considerable attention in the physics 
literature and we believe this duality to be of some interest in its own right. 

A natural starting point for modeling two competing populations is the 
classical Lotka-Volterra model. This is a deterministic model for the evo- 
lution of the total sizes of the two populations, denoted Ni{t), N2{t). They 
are assumed to follow the following system of differential equations: 

where rj, Ki are respectively the intrinsic growth rates and carrying capac- 
ities of the two species and the aij measure the interspecific competition. It 
is easy to check that longterm coexistence of the two populations is possible 
if Ki > ai2K2 and K2 > a2iKi. A number of models have been proposed 
that extend this in two different ways. First, they incorporate spatial struc- 
ture into the populations and second, they assume that the evolution of the 
populations is stochastic. 

It is far from clear how spatial structure affects the chances of longterm 
coexistence for two competing populations. Traditionally, ecologists have be- 
lieved that the local nature of interactions between populations that are dis- 
persed in space promotes coexistence. One reason is the so-called competition- 
colonization trade-off: a weaker competitor that is good at colonization may 
be able to survive by exploiting "gaps" between its competitors. It has also 
been claimed that because in spatial models the population tends to be- 
come segregated into clusters of a single type, the intraspecific competition 
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will be more important than the interspecific competition. Pacala and Levin 
[14] make an attempt to quantify this effect. On the other hand, Neuhauser 
and Pacala [13] propose and analyze a spatial stochastic model for compet- 
ing species in which space actually makes coexistence harder. This suggests, 
then, that in their model, it is actually the interactions at the cluster bound- 
aries that dominate. 

In order to obtain analytic results about spatial stochastic models, sim- 
plifying assumptions must be made. Murrell and Law [12] point out that 
common assumptions are that the parameters of neighborhoods over which 
individuals compete are the same, irrespective of species, or that disper- 
sal and competition neighborhoods are of the same size, but that dropping 
such symmetries can have profound consequences. They argue, using a sim- 
ulation study and the method of moment closure for a specific stochastic 
model in two space dimensions, that spatial structure can promote coexis- 
tence, by showing that in the spatial setting, two populations in which the 
overall strength of interspecific and intraspecific competition is the same 
can coexist, but only if the distance over which individuals sense their het- 
erospecific neighbors (i.e., their competitors) is shorter than that over which 
they sense their conspecific neighbors. They coin the term heteromyopia for 
populations that are "shortsighted" in this way. We explain this concept in 
a little more detail in the context of our first model in Section 2. Although 
this model admits such differences in neighborhood size, our methods are 
not strong enough to confirm the numerical findings of Murrell and Law 
in this context. Indeed, even when the populations migrate in a symmetric 
way and intraspecific and interspecific competition neighborhoods are of the 
same size, although we conjecture (in Section 2) that space does not make 
coexistence harder for our model, our methods are not strong enough to 
provide a rigorous proof of this claim. 

Model I. Following Bolker and Pacala [2], we assume that the strategies 
for survival that individuals in our model can employ are (i) to colonize rel- 
atively unpopulated areas quickly, (ii) to quickly exploit resources in those 
areas and (iii) to tolerate local competition. We take two different popula- 
tions (species) and each can adopt a different combination of strategies for 
survival. In order to simplify the proofs of our results, we suppose our popula- 
tions to be living on Z,'^ (the biologically relevant case is d = 2). The dynamics 
of the model are entirely analogous to those considered by Bolker and Pacala 
[2] and by Murrell and Law [12]. We write {X{t)}t>o = {Xi{t),i G 'L'^}t>o 
and {Y{t)}t>Q = {Yi{t),i € l/}t>a for our two populations. We shall suppose 
that the pair of processes {X{t)}t>o-, \X{t)}t>o solves the following system 
of stochastic differential equations: 

dXi{t) = mij{Xj{t) - Xi{t))dt 
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(2) 



dYi{t)= m[j{Yj{t)-Yi{t))dt 
(3) + a' [w - Y K,Y,{t) - Y ^Ij^jii)] dt 



+ ^aYi{t)dB[{t), 

where {{Bi{t)}t>Q,{B[{t)}t>o,'i S ^''j is a family of independent standard 
Brownian motions. The (bounded nonnegative) parameters ffiij^ ^ij^ ^^j^ 
X[j, 'jij and 7j'j are all supposed to be functions of ||« — i|| alone and to 
vanish for ||i — j|| > R for some R < oo. In other words, the range of both 
migration and interaction for the two populations will be taken to be finite. 

Here, || • || can either denote the lattice distance (so that simple random 
walk is included) or the maximum norm on Z*^, but it will be convenient to 
take the maximum norm. Moreover, notice that by a change of units, there 
is no loss of generality in taking the same a for both populations and indeed, 
hereafter, we may and will set a = 1. 



Remark (Existence and uniqueness). Note that Model I is not covered 
by the now standard results in Shiga and Shimizu [17]. However, Blath, 
Etheridge and Meredith [1] provide, for 1 < p G N, existence of a continuous 
positive (weak) solution in the space = {x £ : ||x||r,4p < oo}, where 
the weighted £P-norm || • ||r,p is defined, for F = {Fjlgd € /^((O, cx))^''), by 
||a^||r,p = (SiGZ<* rj|xj|P)^/*'. We assume Ti/Tj < f{\\i — j\\) for some continu- 
ous function / : [0, oo) [0, oo). For example, set Fj = e""*" for each i S Z'^. 
Uniqueness remains open, after considerable efforts, including those of sev- 
eral experts whom we have consulted. At first sight, one expects to be able 
to prove uniqueness in a suitable weighted ip' space by an application of the 
(infinite-dimensional) Yamada-Watanabe theorem. This works only in the 
special case when Xij and jij both vanish for i^ j. The nonlocal nature of 
the interaction destroys the vestiges of monotonicity available in this special 
case. 



Remark (Blath, Etheridge and Meredith [1]). The full version of this 
paper, Blath, Etherige and Meredith [1], which has also successfully under- 
gone the peer reviewing process of Annals of Applied Probability, proved too 
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long to be published in its entirety. It contains full technical details and 
some additional remarks and is available from the webpages of the authors. 

For the X-population, the first two strategies for survival listed above 
correspond to taking large niij and large aM, while the third corresponds 
to taking small Xij (conspecific competition) and 7jj (interspecific compe- 
tition). By varying M, we can also model how efficiently the species uses 
the available resources: a species that can tolerate lower resource levels will 
have a higher value of M. 

Definition 1.1 (Notions of survival). Let p G [0,1). We shall say that 
the X-population survives for all time with probability greater than p if 
there exists k> such that 

liminfP[Xo(t) > k] >p. 

t— >CXD 

We shall say that both populations persist for all time with probability 
greater than p if there exists k > such that 

¥[yt>0,3i,j £Z'^:X,{t),Yj{t) >K]>p. 

Finally, we shall say that the populations exhibit long-term coexistence with 
probability greater than p if there exists k > such that 

liminfP[Xo(t),yo(i) >k]>P- 

Observe that the third notion is much stronger than the second one. Also, 
note that if 7,^ = 7^^ is zero for all i, j € Z"^, then each population follows 
an independent copy of the so-called stepping stone version of the Bolker- 
Pacala model introduced in Etheridge [9]. There, it is proved that if the 
range of migration is at least as great as the range over which the popu- 
lation interacts with itself (here determined by the {Ajj}), then provided 
that aM is sufficiently large, the population will survive. A partial converse 
of this, proved there only in the context of a continuous-space analogue of 
this model, suggests that this condition is actually necessary, a conclusion 
reached independently by Law, Murrell and Dieckmann [11]. We therefore 
assume from the outset that there exists a constant c > such that for all 
i,j S Z*^, we have mij > cXij (resp. m[j > c\[j) whenever Xij (resp. A-^) is 
strictly positive. Indeed, Theorem 1.5 in Etheridge [9] then tells us that 
if aM > Y^jTUij is sufficiently large (depending on c), the single species 
model for X started from any nontrivial translation invariant initial condi- 
tion survives with positive probability, that is, there exists a k > such that 

liminft_^ooIP [^o(i) > > 0. 

For the competing species model, we will have to make similar and addi- 
tional assumptions. In particular, we shall choose initial conditions in such 
a way that we can find a box where both populations are present, but not 
so prevalent that the competitive interaction between them is too large. 
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Notation and assumptions for Theorem 1.2. 

• The parameters niij, m'^, Xij, X'^j, 'jij and 7^-^ are nonnegative functions 
of II i — j|| alone and vanish for ||z — j || > R for some R < 00. 

• {^ij} {X'^j} are fixed in such a way that there exists a 
constant c> such that, for ah mij,m^j ^ 0, 

(4) -Xij < rriij < cXij and -X'^j < m'^j < cX[y 

For ah i,j such that mij = (resp. m[j = 0) we require Xij = (resp. X[j = 
0). Assume that {ruij} and {m[j} are nondiagonal and of the same range 
and that Xa, A^j > for all i^l/. 

• Let L = max{||j — z|| :mij,m[j ^ 0} < R < 00. 

• Assume that aM > ^j^-^d ruij and a' M' > X^jez^* ''^ij- 

• Let 6 € Z such that max{||j — i\\ ■.'jij or -y'^j 7^ 0} is less than (6 — 1)L. 

• For ?n G N U {00} and < ki < K2 < 00, < n'l < K2 < 00, we write 

(X(0), y(0)) G H{ki,K2; K[,4;m), 

if X{0),Y{0) G 4^ and there exists a box J = {[-m, m]"^ n Z'^} C Z'^, such 
that for all i J, 

(5) Xo{i) e[Ki,K2) and Yo{i) e [k[, K2). 

Remark. One can drop the assumption that the range of {mij} and 
{m'^j} are the same, but this will make the proof much more awkward. 

Theorem 1.2. Under the above assumptions, there exist finite constants 
Mo > 0, Af^ > such that: 

1. for each M > Mq and M' > 0, there is a constant 7 = 'y(M,M') > and 
constants < ki < 00, < K2 < 00, such that if J2j Hj < 7 o.'iT'd 

(X(0), y(0)) G H{ki,oo; 0, 4; (6 + 1/2)L), 

then the X-population survives for all time with probability greater than 
one half; 

2. similarly, for each M' > Mq and M > 0, there is a constant 7' = 7'(M, 
M') > and constants < K2 < 00, < k'^^ < 00, such that if J2j I'ij < l' 
and 

(X(0), y(0)) G F(0, 4, 00; (6 + 1/2)L), 

then the Y -population survives for all time with probability greater than 
one half. 
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Corollary 1.3. Under the conditions of Theorem 1.2, for each pair 
(M, M') with M> Mo and M' > Ml^, there is a pair (7, 7') with 7 > 0, 7' > 



then the X- and Y -populations both persist for all time with positive probabil- 
ity. Moreover, for each such pair, if H[ki, K2] k'^, (6+ 1/2)L) is replaced by 
H{ki,K2; K.i,H2; 00), then there is longterm coexistence with positive proba- 
bility, that is, there exists k > such that 



As we will explain in Section 2, we would conjecture a result very much 
stronger than Theorem 1.2 (or Corollary 1.3). In particular, we provide 
evidence to support the claim that in the biologically relevant case of two 
dimensions, if we take the special case of our model in which a = a' , M = M' 
and lUij = rn[j, then provided that 7^'^- < Ajj (resp. 7jj < A^^) with strict 
inequality whenever Ajj (resp. A^^) 7^ 0, and that the parameters are such 
that if 7jj and ^[j were zero then the single species models would survive, 
then with positive probability, the competing species model will coexist for 
all time. This would be precisely the prediction of the corresponding Lotka- 
Volterra model. If we drop the assumptions a = a' and M = M' , then this 
conjecture must be modified to reflect competition-colonization trade-off. 
We formulate this and other conjectures more carefully in Section 2. In the 
process, we are led to consider our second model of two competing species. 

Model II. Suppose now that the neighborhood over which each individ- 
ual competes is just the site in which it lives, so that the only interaction 
between different points in Z"^ is through migration. In addition, we suppose 
that the migration mechanism for the two populations is the same and that 
the total population size in each site is constant [i.e., Xi{t) + Yi{t) = N > 
for all i gIj'^ and all t > 0]. Let us write Pi{t) = Xi{t)/N for the proportion of 
the total population in i at time t that belongs to the X-population. Then, 
as we will see in Section 2, we arrive at the much simpler model 




liminfP[Xo(t),yo(i) > > 0. 



dpi{t)= ruijiPjit) - pi{t)) dt + .spi{t){l - pi{t)){l - fj,pi{t)) dt 



(6) 




where s 



+ ^N~^Pi{t){l-pi{t))dWi{t) 
aM - a'M' + (a'A-^ - a-iii)N and 



{a'X'ii - a-fii)N + {aXii - a'iii)N 
aM - a'M' + (a'A^j - ayu)N 
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and, finally, {Wi{t),i € Z'^}j>o is a family of independent Brownian motions. 
This model is a system of interacting Fisher- Wright diffusions for gene fre- 
quencies in a spatially structured population. From the results in Shiga and 
Shimizu [17], it follows that if pj(0) £ [0, 1] for all i € Z"^, then this system has 
a continuous, pathwise unique, [0, l]^''-valued strong solution for all times 
t>0. 

If < 1, then in each site i, there is selection in favor of either the X-type 
or the y-type according to whether s > or s < 0. If /U > 1, then in each 
site i, we have selection in favor of heterozygosity if s > and selection in 
favor of homozygosity if s < 0. In the "neutral" case (s = 0), the process has 
a moment dual, the so-called structured coalescent (see, e.g., Shiga [16]), 
and it is easy to show that if d>3, then with positive probability, there 
will be longterm coexistence of our two populations, whereas if d < 2, with 
probability one, eventually only one population will be present. 

Notice that we have selection in favor of heterozygosity precisely when 

{aXii - a'-f'ii)N > aM - a'M' and (a'A-^ - a-fii)N > a'M' - aM. 

We sketch a proof of the following result and present a more detailed analysis 
in a forthcoming work. 

Theorem 1.4. Let {pi{t),i £ Z'^}t>o evolve according to Model II. Sup- 
pose that /X > 1 and let e G (0, 1/4]. Then, ifpi{0) € (e, 1 — e) for some i € Z"^, 
there exists an sq G [0, oo) such that for all s > sq, we have 

P[Vt >0,3ieZ'^:e < pi{t) < 1 - e] > 0. 

Moreover, if Pi{0) G (e, 1 — e) for all f G Z'^, then 

liminf P[e <po{t)<l-e]> 0. 

t— >oo 

In the case when the two populations evolve symmetrically, that is, /x = 2, 
Model II reduces to 

dp^{t) = Y,miJ{p■j{t)-pi{t))dt + spi{t){l-pi{t)){l-2p,{t))dt 

+ ^N-^Pi{t){l-p,{t))dW^{t). 

For general s, there is no convenient dual, but in Lemma 2.1, we find an 
alternative duality with a system of branching annihilating random walks. 

Definition 1.5 (Branching annihilating random walk). The Markov 
process {ni{t),i G Z'^}t>o with values ni{t) G Z_|_ and dynamics described by 

f rij I— fij — 1 . 
< at rate niruij (migration). 
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nj I— > nj + 2 



at rate sni 

at rate ^ni{ni — 1) 



(branching) 



ni<-^ rii — 2 



(annihilation) 



is called a branching annihilating random walk with branching rate s (and 
offspring number two). 

Corollary 1.6. There exists sq>0 such that if s > sq, the branching 
annihilating random walk, started from an even number of particles at time 
zero, will survive for all time with positive probability. 

Remark. Note that in our branching annihilating random walk, a birth 
event results in one individual splitting into three, a net increase of two, 
whereas an annihilation event results in the loss of two particles. As a result, 
we have parity preservation: if we start from an odd number of particles, 
then there will always be an odd number of particles in the system (so, 
in particular, at least one). This is why we restrict the initial number of 
particles in Corollary 1.6 to be even. 

Branching annihilating random walk has received considerable attention 
from physicists (see Tauber [18] for a review). For example, Cardy and 
Tauber [3, 4] consider precisely the process described above. Our conjecture 
for Model II, stated in Section 2, is based on their results, which, in turn, 
are based on perturbation theory and renormalization group calculations. 

The rest of the paper is organized as follows. In Section 2, we explain the 
relationship between our two models and the duality between the symmetric 
form of Model II and branching annihilating random walk. We also make 
some conjectures about the longterm behavior of our two models and relate 
them to results and conjectures for other toy models. The proof of our main 
result will rely upon a comparison with oriented 2A^-dependent percolation 
and so in Section 3, we recall the definition of 2A^-dependent percolation and 
state a suitable comparison result. The proofs of Theorem 1.2 and Corol- 
lary 1.3 are in Section 4 and a sketch of the proof of Theorem 1.4 is in 
Section 4.2.3. Corollary 1.6 will then be an immediate consequence of the 
duality of Model II and branching annihilating random walk. 

2. Heuristics, duality and relation to existing models. 

2.1. Relationship between the two models. Suppose that the evolution of 
our population follows Model I, that is, it is determined by equations (2) 
and (3). We now derive the system of equations governing the propor- 
tion of the total population at time t at site i that belongs to the X- 
subpopulation. We need some notation. If we write Ni(t) = Xi(t) + Yiit) 
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and Pi{t) = Xi{t)/Ni{t), then an application of Ito's formula (and some re- 
arrangement) gives 



dPiit) = E 



N^{t) 



{Pj{t)-Pi{t))dt 



+ E {m,j - m'i^)^^pi{t){l-pj{t)) dt 



(8) 



aM - a'M' + ^ (a'A-^- - a-fij)Nj{t) 



Xpi{t){l-pi{t))dt 

+ ^jNi{ty^pi{t){l-pi{t))dWi{t), 

where {Wi{t),i G Z'^}t>o is a family of independent Brownian motions. 

We concentrate on the case when niij = m[j . Notice that if we also assume 
that Xij,X[j and are zero for j and that the population sizes Ni{t) 

are, in fact, a fixed constant, then we arrive at Model II: 

dPi{t) = Y^viPj('t) - pi{t)) dt + spi{t){l - pi{t)){l - npi{t)) dt 
j 



+ ^N~^pi{t){l-pi{t))dWi{t), 

where s = aM - a'M' + N{a'X'ii - a-fu) and 

{a'X'ii - ajii)N + {aXii - a'-f'ii)N 
aM - a'M' + (a'A^^ - a-/ii)N 



2.2. Conjectures for Model II. Our conjectures for Model II are based 
on the symmetric case, implying that /i = 2. The model then reduces to the 
system (7). In this case, we are able to find a convenient dual process. First, 
we transform the equations. Let Xi[t) = 1 — 2pi{t). Then 



(9) 



dxi{t) = mij{xj{t) —Xi{t))dt 
3 

+ Ux^,{t) - xm dt - ^{i-xmdw,{t). 
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Lemma 2.1. The system (9) is dual to branching annihilating random 
walk with branching rate s/2, denoted {ni(t),i G Z'^}t>o, through the duality 
relationship 

E[x(tp(°)]=E[x(op(*)], 

where x^ = lliezi ■ 

The proof is completely standard (see e.g. Shiga [15]) and is therefore 
omitted. 

Cardy and Tauber [3, 4] consider the branching annihilating random walk 
model of Definition 1.5. In particular, their results suggest that although in 
one dimension the optimal value for sq in Corollary 1.6 is strictly positive, 
in two dimensions one can take sq = 0. 

Conjecture 2.2. For Model II with /i = 2 and d=l, there exists a 
critical value sq > such that the populations described by system (7) will 
both persist for all time with positive probability if and only if s > sq. In 
d = 2, there is positive probability that both populations will persist for all 
time if and only if s > 0. For d>3, this probability is positive if and only if 
s>0. 

Roughly speaking, for d > 2, if there is a homozygous advantage, then 
the population will initially form homogenic clusters, but ultimately it will 
be the interactions at the cluster boundaries that dominate and one type 
will go extinct. In the heterozygous advantage case, there will be long term 
coexistence of species. In one dimension, the heterozygous advantage must 
be "sufficiently strong" if we are to see coexistence. 

In fact, we would go further. In view of the genetic interpretation of Model 
II, it would be odd if the case fi = 2 were pathological, so we expect that in 
d>2, we will have coexistence for any s > 0, fi> 1. 

Conjecture 2.3. Conjecture 2.2 holds true for any fJ, > I, where in 
one dimension, sq will now also depend on /x. 

If this conjecture is true, then in dimensions greater than one, for 

(a'A-j - a-fii)N > a'M' - aM and {aXu - a'jii)N > aM - a'M', 

we have positive probability that both populations survive. Comparing the 
quantities a'X^^ — a'ju and aXu — a'7j'j tells us about the relative effectiveness 
of the X- and y-populations as competitors. If the first is smaller, then the 
X-population is a less effective competitor. However, provided that aM > 
a'M', we can even allow it to be negative and have positive probability 
of survival for the X-population. This reflects a competition-colonization 
trade-off. 
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2.3. Conjectures for Model I. We now turn to Model I. We assume that 
the migration mechanisms governing the two populations are the same. Sup- 
pose first that a = a' , M = M' , Xij = j'^j and X'^j = jij. We then see that the 
system of equations (8) looks like a selectively neutral stepping stone model 
with variable population sizes in each lattice site. If we condition on the tra- 
jectories of those population sizes, then this process will have a dual process: 
a system of coalescing random walks in a space-and-time-varying environ- 
ment. Showing that there is no longterm coexistence of types amounts to 
showing that two independent random walks evolving in this environment 
will, with probability one, eventually meet and coalesce. If the environment 
is sufficiently well behaved, then one might expect this to be true. Problems 
will arise if the environment develops large "holes," so that the walkers never 
meet, or very dense clumps, so that when the walkers do meet, they do so 
in such a heavily populated site that they do not coalesce before moving 
apart again. Much of our proof of Theorem 1.2 is devoted to showing that 
the environment does not clump and a special case of that result says that 
provided both populations are initially present in sufficient numbers in all 
sites, the probability that any given site is in a "hole" at time t is uniformly 
bounded below. We therefore conjecture that in the neutral case. Model I 
will behave qualitatively in the same way as Model II. In the biologically 
relevant case of two spatial dimensions, we have been unable to produce a 
proof. 

More generally, we believe, still assuming that niij = m'^j and a = a' , 
M = M' , provided that at least one population persists, the question of 
longterm coexistence of the populations described by Model I will not be 
changed by assuming that competition only acts within individual lattice 
sites and, moreover, in that case, the question of coexistence will be the 
same as for the populations described by Model II. Namely, we make the 
following conjecture. 

Conjecture 2.4. Let niij = m[p a = a' , M = M' he fixed. Suppose that 
both X- and Y -populations start from nontrivial translation invariant initial 
conditions and that the parameters are such that each population has positive 
chance of survival in the absence of the other. Further, let Xij = X[j,^ij = ^[y 

1. If Xij < jij for all j , then eventually, only one population will be present. 

2. // Xij > 'jij for all j , then if d> 2, with positive probability, the popula- 
tions will exhibit longterm coexistence. In one dimension, the same result 
will hold true provided that Xij — jij is sufficiently large. 

3. If Xij = ^ij and d>3, then with positive probability, both populations will 
exhibit longterm coexistence. 

If d<2, then with probability one, one of the populations will eventu- 
ally die out. 
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When aM ^ a' M' , we would once again expect to see a competition- 
colonization trade-off. 

2.4. Heteromyopia. In view of equation (8), it is easy to see that Mur- 
rell and Law's heteromyopia might lead to coexistence. They work in a 
continuous space with the strength of competition between individuals de- 
caying with their distance apart, according to a Gaussian kernel. The ana- 
logue of their model in our setting is the symmetric version of Model I with 
\j = A(||i — ill), 7ij =7(11* — ill), where the functions A and 7 are mono- 
tone decreasing and J2j ^ij = J2j lij ■• but the range of \ij is greater than 
that of 7ij. We can think of the effect of this as follows. Over small scales, 
we have homozygous advantage, over larger scales, heterozygous advantage. 
Again, we expect to see the population forming homogenic clusters, but now 
the cluster boundaries will be maintained because the heterogeneity there 
confers an advantage to individuals within the clusters which counteracts 
the disadvantage to the individuals actually on the boundary. Reversing the 
sign to give populations with "heterohyperopia" produces the opposite ef- 
fect. This is not stable and Murrell and Law observe founder control in this 
case, which means that the outcome of the competition is entirely deter- 
mined by the initial conditions. 

Conjecture 2.5. Assume nontrivial translation invariant initial con- 
ditions. Suppose that the parameters of Model I are symmetric and such that 
in the absence of the competitor, each population survives for all time with 
positive probability. For d>2, if the populations are heteromyopic, then we 
will see longterm coexistence with positive probability, whereas in d=\, the 
populations must be strongly heteromyopic for there to be longterm coexis- 
tence with positive probability. 

2.5. Relation to existing models. 

The Murrell-Law model. Our conjectures for Model I are entirely in 
agreement with the numerical results of Murrell and Law [12]. They analyze 
a stochastic version of a continuous-space Lotka-Volterra system, similar to 
ours. The evolution is characterized in terms of moment equations. These 
moment equations were derived from a stochastic individual-based model 
by Dieckmann and Law [7] . Although the assumption of a spatially contin- 
uous environment is clearly desirable, the price that they pay is that there 
are very few analytic tools available for the study of the resulting popula- 
tion models, so they use moment closure, assuming in this case a "power-1" 
closure. In particular, they ignore dynamics of all spatial moments beyond 
order two. In view of the clustering behavior that is characteristic of pop- 
ulations evolving according to spatial branching models in two dimensions. 
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this method has potential pitfalls. In fact, the control of the clumping of 
the populations that forms an essential part of our proof of Theorem 1.2 
also adds considerable credibility to the moment closure technique for these 
models and hence to the numerical predictions of Murrell and Law. 

Neuhauser and Pacala [2] also consider an explicitly spatial stochastic 
version of the Lotka-Volterra model. Their model is lattice based, but, in 
contrast to ours, allows only a single individual to live at each lattice site. 
Moreover, there is instant recolonization so that there will always be exactly 
one individual at each site in Z"^. This fixed population size makes it more 
analogous to Model II than to Model I. 

Definition 2.6 (Neuhauser-Pacala model). The Markov process {r]i{t), 
i G Z,'^}t>o in which 7]i{t) € {1,2} and with dynamics: 

1. if r]i{t) = 1, it becomes 2 at rate xfj+fi (/i + '^12/2)5 

2. if r]i{t) = 2, it becomes 1 at rate xf2+fi (-^^ ~'~ '-'^21/1)5 

where fk{i) = '^^^ "''^^j^^'f and Mi = i + {j:0< < R}, will be said 

to follow the Neuhauser-Pacala (stochastic spatial Lotka-Volterra) model. 

The idea is that an individual of type k will die at a rate determined 
by the proportion of its neighbors that are conspecific plus some constant 
multiple of the proportion of heterospecific neighbors. Thus, for example, if, 
in Model I, we took Xij and 'jij to have the same range and to be constant on 
that range, then a small value of ai2 would correspond to the ratio 'Jij/Kj 
being small. The dead individual is immediately replaced by an offspring of 
one of its neighbors chosen according to a weight that reflects the relative 
fecundity of the two types. Thus, for example, A > 1 would reflect type 2 
being more fecund than type 1. In Model I, this would be modeled by taking 
a'M' > aM. Let us recall some results for this model. 

Theorem 2.7 (Neuhauser and Pacala [13], Theorem 1). Suppose that 
X= 1, d = 1 or 2 and 012 = 021 = ct- 

1. If a = 0, then, except for the one- dimensional nearest-neighbor case, prod- 
uct measure with density 1/2 is the limiting distribution starting from any 
nontrivial initial distribution. 

2. If a is sufficiently small (depending on R), then coexistence is possible 
except for the one- dimensional nearest-neighbor case. 

For Model I, a result entirely analogous to part (2) is a special case of 
Theorem 1.2. For Model II, the analogue is Theorem 1.4. If we believe Con- 
jecture 2.2, then although in d = 1 we require the condition "a sufficiently 
small," in d = 2, the corresponding result is true for all a <1. This corre- 
sponds to Conjecture 1 of Neuhauser and Pacala [13]. 
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Fig. 1. The persistence region of the Neuhauser-Pacala model. 



Theorem 2.8 (Neuhauser and Pacala [13], Corollary 1). Suppose that 
A = 1. Write n = \J\f\ for the number of lattice sites in a neighborhood. 
Species 1 competitively excludes species 2 if 



In particular, this result shows that the values of (0:12,021) ^or which 
both populations persist for all time are contained in the shaded region 
in Figure 1. This is a reduction from the range of values predicted by the 
mean field model. The case A = 1 corresponds, in our setting, to taking 
iTLij = m'^j, a = a' and M = M' , so in view of Conjecture 2.3, we expect 
that the coexistence region for Model II in two dimensions corresponds to 
the whole region [0,1) x [0,1) in the (012, a2i)-plane, that is, the region 
predicted by the mean field model. 

Cox and Perkins [5] show that a sequence of processes following the 
Neuhauser-Pacala model, when suitably rescaled in space and time, con- 
verges to a superBrownian motion with a nontrivial drift. In low dimensions, 
they restrict to long-range models, whereas in dimensions d>3, they can 
also consider the nearest-neighbor case. SuperBrownian motion has emerged 
as a universal limit of critical spatial systems above the critical dimension 
and these results can be seen as special cases of a general convergence the- 
orem for perturbations of the voter model. In a recent preprint. Cox and 
Perkins [6] show that in dimensions d>3, the drift in the superBrownian 
motion is connected to questions of coexistence in the Neuhauser-Pacala 
model. Using this connection, they obtain additional information about the 




for 021 G 1 , 1 , 

V n 

for 021 > 1- 



Species 2 competitively excludes species 



1 ^/ 




for 021 € (0,1] 
for 021 > 1. 
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parameter regions in which survival of one type (resp. coexistence) holds. 
The biologically relevant case d = 2 is a topic of their current research. 

3. 2Ar-dependent oriented percolation. We now turn to proving our re- 
sults. Since our proofs will rely upon comparison with 2iV-dependent ori- 
ented percolation, we first briefly recall some well-known facts which can 
be found, for example, in Durrett [8]. The insistence on 2N- instead of A^- 
dependent percolation will be explained in the remark below Theorem 3.5. 

Oriented percolation will be defined on the lattice 

C := {{x, n) G : x + n is even, n > 0}. 

This set is made into a graph by inserting edges from (x, n) to (x + 1, n + 1) 
and to (x — l,n + 1). It is convenient to think of n as time. We introduce 
a family of {0, l}-valued random variables uj{x,n) at sites (x,n) € C A site 
(x,n) is called open if uj{x,n) = 1 and closed if a;(x,n) = 0. Given such a 
family of random variables and integers < m < n, we say that (y, n) £ C can 
be reached from {x, m) if there is a sequence of points x = Xm, Xm+i, ■ ■ ■ ,Xn = 
y such that \xk — Xk~.i \ = 1 and w(xfc, k) = 1 ioi m < k < n. We write this as 
(x, m) — > (y, n). Finally, given an initial condition Wo C 2Z = {x : (x, 0) € £}, 
we may define a percolation process {Wn}n>o by setting, for each n > 0, 
Wn = {y '■ {x, 0) — > (y, n) for some x G Wo}. 

Definition 3.1. Let 9 G (0, 1) and N eN. We say that an oriented per- 
colation process {Wn}n.>05 determined by ra)}(a;,n)G£5 is 2N -dependent 
with density at least 1 — ^ if, for any finite set of indices / such that | (x^ , n^) 
{xi,ni)\\ >2N for all k^l€ I, we have P[a;(xfc, n^) = 0, G /] <6 



i\ 



Define Cq = {{y, n) G : (0, 0) — > {y, n)} as the open cluster containing the 
origin. We say that percolation occurs if |Co| = oo. We first cite a result which 
gives us a lower bound for the probability of percolation depending on 9 and 
N. A proof can be found in Durrett [8]. 

Theorem 3.2. If 9 < q-H^n+i^ ^ ^/^g^ py^-^j < oo] < 559^^^^+^^)^ < ^. 

For particle system models of evolving populations, a standard strategy 
for showing survival is to construct a suitable coupling with oriented per- 
colation. Our approach amounts to a modification of this strategy to cope 
with the interactions between the two populations, so we now describe the 
relevant comparison theorems. Once again, we are citing Durrett [8], but 
we also present a modified version of the results which are adequate for our 
purposes. In Durrett's terminology, let {^^(n), z G Z'^}n>o denote a trans- 
lation invariant time-homogeneous finite-range flip process with state space 
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17 = {0, 1} , constructed from the usual graphical representation. Let L e N 
be fixed. We write 

(10) H={^e {0, 1}^' : = 1 Vi G [-L/2, L/2f n Z"^} 

and m • H for the translation of H by some integer m with respect to the 
first component, that is, 

(11) m-H = {Ce{0, 1}^' : = 1 Vi G mLei + [-L/2, L/2]'^ n Z'^}, 
where ei is the unit vector in the direction of the first component. 

Definition 3.3. Fix N eN and 9 € (0, 1). We shall say that the pro- 
cess {C{n)}n>o fulfils the classical comparison assumptions for and 9 
if, for each configuration ^ G ff, there exists a "good event" G^, measur- 
able with respect to the graphical representation of the flip process inside 
[-NL,NL]'^ X [0, 1] and with F[G^] > 1 - 6*, so that if ^(0) = ^, then on G^, 
we have ^(1) G (+1) • H n (-1) • H. 

Unfortunately, the flip processes we are going to consider in the next sec- 
tion are functionals of more general underlying stochastic processes driven 
by independent Brownian motions and cannot be obtained from a graphi- 
cal representation. However, a modified comparison result holds, based on 
the behavior of the Brownian increments. More explicitly, we are going to 
construct our flip process in terms of the system of stochastic differential 
equations X,Y from Model I. Define {^i{n),i G Z'^}n^^ by 

{1, if Xi{2n)>ci and Yj{2n)<C2 
Vj Gi + [-6L,6L]^nZ'^, 
0, otherwise, 

where ci , C2 are finite positive numbers and b and L satisfy the assumption 
of Theorem 1.2. Note the different time scales for ^ and X,Y (the usefulness 
of this time change will become clear in the next section) and observe that 
X, Y are time-homogeneous and the underlying system of driving Brownian 
motions is translation invariant. 

Definition 3.4. Assume N ■.= {b + 2) and 9 e {0,1). Define the events 
H and m ■ H for some integer m €Z in terms of the process {(,{n)}n>o in 
the same way as in (10) and (11). Define 

(13) J^*iNL, [0, 2]) = a{Bi{s),B'^it) : s, t G [0, 2];i,j G [-NL, NLf n Z^}, 

where {{Bi{s)}s>o, {Bj{t)}t>o,i,j G Z'^} is the family of independent stan- 
dard Brownian motions as in Model I; see (2) and (3). We shall say that 
the process {'^(?^)}n>o fulfils the modified comparison assumptions for NL 
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and 9 if, for each configuration ^ E if, there exists a "good event" G 
T*{NL, [0, 2]), with F[G^] > 1 - 6*, such that if ^(0) = ^, then on G^, 

(14) ^-(i)e(+i).ijn(-l)-i7. 

In other words, if ^ has all I's in the box of side length L about the 
origin at time 0, then at time 1 (measured in time units for the ^ process), 
with probability at least 1 — 0, it has successfully "invaded" the boxes of 
side length L translated by —Lei and Lei in a way that is measurable with 
respect to the Brownian increments inside the box around the origin of side 
2NL and up to time [0,2] (measured in time units for X,Y). 

The following two classical theorems (see, e.g., Durrett [8]) apply in both 
of the above settings and complete this section. 

Theorem 3.5. If the classical (resp. modified^ comparison assumptions 
hold for ^ (resp. ^) for some S N and 9 € (0, 1), we may define random 
variables uj{x,n) such that Xn := {(m,n) E C:^{n) G m ■ H} [resp. Xn ■ = 
{{m,n) E C:^{n) E m ■ H}] dominates a 2N -dependent oriented percolation 
process {Wn}, defined on C, with initial configuration Wq = Xq (resp. Xq) 
and density parameter at least 1 — 0, that is, Wn C Xn (resp. Xn) for all 
n E N, where Wn = {y ■ {x, 0) — > {y, n) for some x E Wo}. 

Theorem 3.6. Suppose that {Wn}n>o is a 2N -dependent oriented per- 
colation process, started from the trivial initial state Wo (2;) = 1 for all x. If 
0< 6-4(4^+1)', then 

liminfP[OE W2n] > i- 

Eventually, this theorem will be the key to proving the coexistence (with 
positive probability) of {X,Y) in Model I and Model II. 

4. Proofs. The key to the proof of Theorem 1.2 is to consider the Y- 
population as providing a random environment in which the X-population 
evolves. Of course, the environment itself depends on the X-population, 
but we obtain some control of the behavior of the environment that is in- 
dependent of the evolution of the X-population. This "decoupling" (and a 
symmetric argument for Y) then reduces the coexistence problem to that 
of survival of a single population: if the X- and y-populations can each 
be shown to survive for all times with probability greater than one half, 
then longterm coexistence (with positive probability) will follow. We attack 
the question of survival of the X-population (resp. y-population) by com- 
parison with a 2A'^-dependent oriented percolation process. To this end, we 
establish the existence of the corresponding "good events" as required in 
Definition 3.4. 
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4.1. A spin system and estimation of related flip probabilities. The main 
step is to construct two spin systems, one for each of the X- and y-populations, 
that play the role of {^i{n),i G Z'^}„>o of the last section for some suitable 
constants ci, C2. Indeed, we consider the spin system {Q{n) , rii{n) , i € Z'^}„>o, 
where {C,i{n),i G Z'^}„>o is defined by 



(15) Q{n) -- 
and similarly, 

(16) r]i{n) 



M 

1, if Xi{2n) > — and Yj{2n) < a'M' 

yj £i + [-bL,bLYnZ'^, 
0, otherwise. 



M' 

1, if Yi{2n) > — and Xj{2n) < aM 

ei + [-bL,bLfnZ'^, 
0, otherwise, 



where K := 2aMc + 1, K' := 2a'M'c + 1 and a, a' are finite positive con- 
stants to be determined later [see (44) in the proof of Lemma 4.5 (resp. the 
symmetric result for a)]. Recall that L denotes the range of the intraspecific 
interaction and b denotes the smallest positive integer such that the range 
of {7jj} (resp. {"flj}) is less than (6 — 1) • L. 

With these definitions, one expects that if the system {Cin'),r]{n)}n>o ex- 
hibits longterm coexistence in discrete time, then the system {X{t),Y{t)}t>o 
exhibits coexistence in continuous time and, in fact, this will follow from our 
proof. The convenience of the time change ni-^ 2n in (15) and (16) will be- 
come clear when carrying out the comparison arguments in Section 4.2.1. 

Outline of this subsection. The notation from Section 4.1.1 is necessary to 
define the suitably measurable "good events" G(^,Gr^. Section 4.1.2 provides 
"flip probabilities" related to the spin system (" via comparisons in terms 
of the aforementioned behavior of the one-dimensional diffusions, under the 
additional condition that the system X evolves in a "safe environment," that 
is, given some bounds on the local Y-population. Finally, in Section 4.1.3, 
we will find conditions so that the "safe environment" assumption holds 
for the y-population in a way that is independent of the evolution of the 
X-population, again by making use of comparisons to one-dimensional dif- 
fusions. 

4.1.1. Some notation and lattices of one- dimensional diffusions. 

Definition 4.1. Let {{Bi{s)}s>o,{Bj{t)}t>o,i,j&'^'^} be the family of 
independent standard Brownian motions driving Model I. Fix i G Z''. For 
n G N, n > 0, define 

!F{i, n, u) := a{Bi{n + s) — Bi{n) : s G [0, u]} 
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and define J-''{i,n,u) accordingly in terms of B' . Moreover, let 
T{i,NL,n,u) ■.= a{Bj{n + s)-Bj{n):se [0,u], j e i + [-NL,NLf nZ'^}, 
and similarly define J-' {i,N L,n,u) in terms of B' . 

These a-algebras will be used to construct suitably measurable events 
G^,Gfj. Recall from Definition 3.4 that 

J^* {NL, [0,2]) =T{0,NL, 0,2) V T'{0,NL, 0,2). 

Definition 4.2. Let i gZ'^ and assume that the constants a', M,\,JJ > 
are chosen such that XU > 2M. Moreover, let Di,D2 > 0. Then, for each 
J € i + [—NL,NLY n Z'^, define the one-dimensional diffusions {Zj{t)}t>o, 
{Zj{t)}t>o, {^ji't)}t>o and {Zj{t)}t>o, driven by independent standard Brow- 
nian motions {Wj{t)}t>o, by 

(17) dZj{t) = a'{M - XZj{t))Zj{t) dt + ^Zj{t)dWj{t) 
(logistic Feller diffusion), 

(18) dZj{t) = DiZj{t) dt + ^JZj{t)dWj{t) 
(supercritical Feller diffusion), 

(19) dZj{t) = D2dt + DiZj{t) dt + ^JZj{t) dWj{t) 
(supercritical Feller diffusion with constant positive immigration), 

(20) dZj{t) = a'{M - \U)Zj{t) dt + ^Zj{t) dWj(t) 
(subcritical Feller diffusion). 

Since each of the four diffusions admits a (continuous) unique strong solu- 
tion, we may assume them to be driven by some given family of independent 
Brownian motions, in particular, those obtained either from (2) or from (3). 

4.1.2. Injection and recovery probabilities for the X -population. Suppose 
that we are interested in the behavior of the X-population within the time 
interval [n,n + 1] at site i and that we already know [recall (15)] that 

(21) max Yi(t)<2a'M' \/j (£i + [-bL,bLf nZ'^. 

te[n,n+l] 

Assume the interspecific competition {'Jij} is chosen such that (21) implies 

(22) max V 7^7 !!(*)< 1 \/j £ i + [-L/2, L/2f nZ'^. 

t£\n,n+l\ 
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This is possible since the range of {7^7} is, by assumption, less than (6 — 
1)L [choose, e.g., Y.j1ij < (2a'M')"^]. We will later (in Section 4.1.3) con- 
struct events that are measurable with respect to either J^'{i,NL,n,2) or 
{i, NL,n — 1,2), which imply (21) and are of sufficiently high probability 
[cf. (39) and (40)]. For the moment, to aid intuition and to simplify notation, 
we will say that, in either case, a suitably measurable "sa/e environment con- 
dition G'^^^{i,n) holds at site i and time re," which implies that (21) [and for 
sufficiently small 7ij, also (22)] holds and which will be explicitly determined 
later. 

We now consider "flip probabilities" for the X-population that are closely 
linked to the flip probabilities of the ^-population, introduced in (15) and 
(16), under the "safe environment condition" Gsec(^)^) sif^ i and at time 
re e N. 

Lemma 4.3 (Infection and nonrecovery) . Let re € N and i ^TJ^. Let 
a, {reijj}, he fixed. Given the parameters for the Y -population and some 
a' > 0, choose {^Jij} such that (21) implies (22). Then, for any e E (0,1), 
there exists a finite constant Mq > such that if M > Mq and K : = 2aMc-\- 
1, for each j G i + [— L/2, L/2]°' n Z'^, there exist events 

(23) GnonTcc{i,n) € J^{i,n,l) and Ginfcc(i, J, re) G :r(i, L/2, re, 1), 

both measurable w.r.t. T{i,NL,n,l), such that the following holds: 

(i) we have 

{{Xi{n) > M/K} n G^ec(^,^) n Gnonrec(i,n)} C {Xi{n + 1) > M/K} 
( ''nonrecovery'' ) and for the "nonrecovery probability" 

Pnonrec 

{i,n), we have 

the bound 

(24) Pnonrec(«,n) := P[Gnonrec(i, ■«)] > I - E] 

(ii) moreover, 

{{Xi{n) < M/K} n {3j:mij > 0,Xj{n) > M/K} n G',^^{i,n)} 
C{Xi{n+l)> M/K} 
(''infection by an occupied neighbor" ) and for the "infection probability" 

Pinfec 

{i,j,n), we have the bound 

(25) Pmiec{hj,n) :=P[Ginfec(^,J»] > 1 - £• 



Proof, (i) We distinguish the two cases Xi{n) G {M/K,{3/2)M/K) 
and Xi{n) > (3/2) M/K. 
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Case 1. Suppose that Xi{n) > {'i/2)M / K and introduce the first hitting 
time of level M/K from above after time n: 

(26) ^M/Ki^) ■■= inf{i >n:Xi{t) = M/K}. 

Our goal is to establish the existence of a suitably measurable event Gnonrcc(^) n) E 
J^{i,n,l) so that Gnonrccii,n) implies, under the above conditions, that 
r^y^(n) > 1. To this end, we set up a suitable comparison to a one-dimensional 
diffusion. Indeed, rearranging the drift in equation (2), as long as Yj{t) < 
a'M' for all j ei + [-bL, bVf n Z'^ and hence 

Y^lMt) < 1 for alii ei + [-L/2,L/2]'^nZ'^ 
I 



holds, and as long as Xi < 2M/K, we have 

~K 



^-^ / 2M \ 

dXi{t) > ruij - a—Xij )Xj{t) dt 



(27) 

+ laM - mij - a] X^{t) dt + ^/x~(t) dBi{t) . 

We now check that the first component of the drift on the right-hand side 
is positive. Indeed, from the assumption (4), we obtain 

2M ^ 2M 
(28) rriij - a——Xij > rriij - a——c'mij, 

K K 

which is positive by our choice of K = 2aMc + 1, for all j G l/. Moreover, we 
have, for each M > 1, that M/K = M/{2aMc+l) e ((2ac-M)-^ (2qc)-^). 
Under these conditions, (27) implies that 



(29) dXi{t)>\aM - irnj - a\Xi{t) dt + ^J Xi{t) dBi (t) 

and so while Xi S [0,2M/K], using Corollary A. 2 to the Ikeda-Watanabe 
Comparison Theorem A.l (both found in the Appendix), we may compare 
Xi to a dominated supercritical Feller diffusion Zi defined in (18), with 
initial value Zi{n) := {3/2)M/K and Di = Di{M) = {aM - Y.jezd ^ij - ")> 
and driven by the same Brownian motion, that is, {Wi{t)}t>o ■= {Bi{t)}t>o- 
It is an important observation in Corollary A. 2 that actually more is true: 
our domination argument holds not only up to the time when Xi leaves the 
interval [0,2A'I/K] for the first time, but, in fact, for as long as Zi takes 

values inside this interval, that is, up to the first exit time t^Ij^j^ [defined 
as in (26)]. 
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Note that for AI > the "supercriticality" (i.e., positive drift) 

Di = (aM — ^niij — a) in (29) tends to oo as M — > oo, while maintaining 
the condition Xi{n) > l/(2ac). 

We now make use of the comparison. Indeed, for t>n, as long as Zi{t) 
stays inside the interval [0,2M/K] and given that initially Xi{n) > Zi{n) := 
(3/2)M/ir, we have that Xi dominates Zi. To obtain a comparison that is 
valid throughout the whole time interval [n,n + 1], we go one step further 
and modify Zi so that whenever Zj hits level 2M/K (and thus is about to 
leave the area in which the comparison holds true) , we restart the process Zi 
at level {'i/2)M / K and repeat this procedure as often as necessary, so that 
the comparison holds for all times i € [n,n + 1]. More precisely, we define a 

sequence of stopping times, beginning with i^2M/k^^^ ■~ '^2M/k^''^^^ restart 

the Zi process at this time, setting Zj(z^^'j-y^(n, 1)) := {3/2)M/K , and then 
iterate this procedure, considering, for m G N, 

(30) '^2AI/K m+l):= mi{t > vll^^^ (n, ?n) : Z^ (t) = {2M)/K}, 
and again restarting the Zi process accordingly, that is, setting 

Uu^M/Ki^,m+l)) = (3/2)(M//0. 

Note that ^'^m/k^'^-' "*-) T oo a.s. as m ^ oo. For definiteness, set t^2M/k^'^-' 
n. We define the i.i.d. positive lengths of the corresponding upcrossing in- 
tervals from (3/2) M/K to 2M/K for m > 1 by 

(31) fm ■■= U^M/K^'^^^) - ^2M/k('^^^- l)' 

Now observe that there is an event ^^^^.^(.(i, n), defined only in terms of 
the Brownian increments {Bi{n + s) — Biin) :s £ [0, 1]}, and hence being an 
element of J-{i, re, 1) , such that if we start our modified diffusion Zi in Zi (n) = 

{3/2)M/K, this event Gl^„^j.^^{i,n) actually equals {''"^/^(re) > re + 1}. [The 
set Gnojn,(5(,(i, re) contains all such uj, such that the corresponding Brownian 
increments lead to the desired behavior if they drive the modified diffusion 
Zi started at time re in (3/2)M/K.] Moreover, by our comparison, the event 

{rj^'y^(n) > re + 1} implies that {Tj^/^(re) > re + 1}, which, in turn, implies 
that Xi{n + 1)> M/K. 

It remains to show that the event Gl^^^j.^^.{i,n) has sufficiently high prob- 
ability. To this end, note that the number of upcrossings of the modified and 
suitably restarted process Zi from level {3/2)M / K to level 2M/K before 
the first downcrossing from {3/2)M / K to M/K is a geometric random vari- 
able with positive parameter qM,K '■—^['''m/k^''^) ^ standard 
speed measure and scale function computations, it is possible to show, for 
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sufficiently large Di and hence M, that qM,K ^ exp{—DiM/K); see Blath, 
Meredith and Etheridge [1] for details. Moreover, by a similar computation, 
the expected time for Zi to exit from the interval {M/K,2M/K) when be- 
ing started in {?,/2)M/K is bounded below by for sufficiently large Di. 
Hence, for such Di, 

(32) nTi]>nrk,iK.2MlK^> ^ 



'{M/K,2M/K)\ - ^2Di' 

Now let T> denote the number of upcrossings before the first "success," that 
is, a downcrossing from {2>/2)M / K to M/K. For each € N, we may then 
write 

P[r5/^(n) < n + 1] = P[r|;/^(n) <n + l-V<N] 



< ¥[!) <N] + 



V 

Y,fi<l;V>N 

li=l 



- N 

<Nexp{-DiM/K) + F Y,Ti<l 

.1=1 

using Bernoulli's inequality. Since, by (32), for large Di, the expectation of 
the length of the i.i.d. upcrossing intervals {T,} of the modified and suit- 
ably restarted process is bounded below by l/{32Di), the number of such 
upcrossing intervals up to time 1 is at most of order Di. Hence, by the Law 
of Large Numbers, we can find a constant d such that for N : = d- Di and all 
sufficiently large Di, the last term on the right-hand side is bounded by e/4. 
Since the first term on the right-hand side still decreases exponentially in 
Di once d is fixed (the linearly increasing prefactor being squashed), for Di 
and hence M sufficiently large, this bound holds simultaneously for the first 
and the last term and we arrive at the desired_result: under the above con- 
ditions, with Zi{n) = {3/2)M/K, we have P[rf//^(n) <n + l]<e/2 for Di 

and hence M sufficiently large, which in turn implies F[Gl^g^j.^^{i, n)] > 1 — |, 
so that Case 1 of part (i) follows. 

Case 2. Now suppose that M/K < Xi{n) < {3/2)M/K. In this case, 
we cannot find a uniform lower bound on the probability on the previ- 
ously considered event {T^Ji^in) > n + 1}, and hence on the probability of 
{rj^'^^(n) > n + 1}, that is sufficiently large. However, we may still use the 

same comparison as above to a dominated supercritical Feller diffusion Zj, 
so that the comparison works as long as Zi stays below 2M/K. This time, 
set Zi{n) = M/K < Xiin) and observe that there is a constant Mq > such 
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that for all M > Mq, the deterministic drift in the supercritical Feller diffu- 
sion Zi will achieve two goals with sufficiently high probability: first, make 
Zi hit level {3/2)M/K within the time interval [n,n + 1/2] with sufficiently 
high probability and second, after hitting level {3/2)M/K , arguing just as 
in the first part of the lemma, ensure that there will be no further down- 
crossing from {3/2)M/K to M/K up to time [n + 1]. Thus, once again, we 
can find a measurable event G'^^^^^^{i^n) S J^{i,n, 1), depending only on the 
corresponding Brownian increments, so that given Zi{n) = MjK^ by com- 
parison, Gnonrec(^5 implies that Xiin ^ 1) > M/K and, moreover, that 
I^[^nonrec(*! ^)] > 1 ~ Hence, the result also holds in Case 2. Finally, in 
view of both cases, choose 

G'nonrcc(«, n) := Gl^^j.^^{i, n) n Gnoni.cc(^> ^ -^(^i '^^ 1) ^ NL, n, 1) 

and part (i) follows. 

To prove part (ii), we begin with some preliminary considerations. Note 
that by using the same comparison and similar arguments as before, again 
considering suitable up- and downer ossings, this time from M/K down 
to M/{2K), we can actually go one step further and find a finite con- 
stant Mq such that if M > Mq, and for j e such that rriij > 0, there 
exists an event Gpcr-occ{hj,n') £ J^{i, NL,n,l), such that given Xj{n) > 
M/K and Gg^^.{i,n), the event Gpcr-occ{hj,n) implies Xj{n + 1) > M/K 

and 'Tjy//(2K)(''^^ > n + I and, moreover, we have 



Note that, once again, we use the assumption that the range of the {'Jij} 
is less than (6 — 1)L so that the "safe environment condition," in particular 
(22), allows comparisons of the above type also for site j. 

We are now prepared to consider the infection probability at a site i in the 
presence of at least one occupied neighbor, say, at j* . Again, assuming (4), 
we use a comparison based on Corollary A. 2. This theorem can be found in 
the Appendix. This time, we rearrange the drift so as to highlight the role of 
immigration of mass to an unoccupied site from occupied neighbors. Once 
immigrated, we can then compare the evolution of the mass to a supercrit- 
ical continuous-state branching process, as before. Indeed, considering the 
drift in equation (2), observe that as long as Xi{t) < 2M/K and given the 
existence of at least one neighbor at some site j* G Z,'^ with rriij* > and 
Xj*(n) > M/K (noting that rriij* is bounded below by some d > since the 
family {niij} is of finite range), we have that, as long as t satisfies 



(33) 



■occ 



(^,j»]>l-e/2. 



(34) 



n<t <t' 




in), 
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by our choice K = 2aMc + 1, 

dxm > ~ ^^'"^^"'^ dt 

4ac + 2 

(35) 

+ iaM - J2 rn^j - a] Xi{t) dt + y/x~{t) dBi{t), 

assuming M > max{l, ^j-mij/a}. Thus, at the uninfected site i, after time 
n, as long as (34) holds, we may compare the evolution of the process Xi 
to a dominated supercritical branching process Zi with constant immigra- 
tion, as defined in (19), and driven by the same Brownian motion, that is, 
{Wi{t)}t>o := {B^{t)}t>o, where D2 = {rriij* - Aij./c)(4ac + 2)-^ > and 
Di = {aM — J2j i^ij — a) > 0. Note that we here use the fact that we have 
assumed the strict inequality ^Ajj* < rriij* from (4) to obtain a strictly pos- 
itive immigration rate. Also, note that the rate of immigration is bounded 
below independently of M. Again, the "supercriticality" {aM — J2j "f^ij ~ oi) 
tends to 00 as M tends to 00. Hence, arguing as before in part (i), this time 
starting the dominated process in Zi{n) = < Xi{n), stopping Zj once it 
reaches level 2M/K and then restarting at {'i/2)M / K if necessary, we may 
find a constant Mq > such that if M > Mq, under the above conditions, 
the event 

Zi{n + 1)> M/K imphes Xi{n + l)> M/K 

and has probability greater than e/2. As before, there is an event G^^^^^{i,j*,n) 
!F{i,n,l), depending solely on the Brownian increments at i within the 
time interval [n,n + 1], so that if we start our modified diffusion at time 
n in Zi{n) = 0, driven by the corresponding Brownian increments, we have 
GUeci^d^n) = {Zi{n + 1) > M/K}. Moreover, 

P[Gr„fc,(i,r,n)]>l-e/2. 

Now, observe that, due to our preparations, (34) is actually being guaranteed 

by Gper-occ(^) 

j*,n) up to time n + l. Combining both events (which are 
actually independent), we define 

Ginicc{i,j*,n) := G*^{^^{iJ*,n) n Gpcr-occ{i,j* ,n) e L/2,n,l) 

and finally see that P[Ginfec(^i J*) "■)] > 1 — £• Together with the fact that 
given Xj*{n) > M/K, our event Ginfcc(^i j)"') implies Xi{n -|- 1) > M/K, 
defining Mq = max{MQ , Mq , Mg , Mq } completes the proof of Lemma 4.3. 
□ 
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4.1.3. Control of the environment. We now find estimates for suitably 
measurable events, a combination of which wih later provide the "safe en- 
vironment condition" Gg^^{i,n) at some site i G Z"^ and time N, which 
wiU, in turn, imply (21). Again, this is done via suitable comparisons to 
one-dimensional diffusions. 

Lemma 4.4 (Control of the environment). Let n S N and i S Z*^. Let 
a' , {m'^j} , {Xij} be fixed. Then, for any e S (0,1), there is a finite constant 
Vq>0 such that for all v' > v'q, there exists an event E'{i,n)[v'] S J^'{i,n,2) 
such that for all a'M' > J2j ^'ij! 

F[E'{i,n)[v']] > 1-e. 

Moreover, 

(36) E'{i,n)[v'] C | sup Yi{n + 1 + s) < 2v'M'\ 

lo<s<l J 

and 

(37) {E'{i,n)[v']n{Y,{n)e[0,v'M']}}cl sup Yi{n + s) < 2v' M'\ . 

lo<s<2 J 

These results hold true for any choice of {7j-j} and any parameter values for 
the X -population in Model I. 

Remarks. 1. In the proof, we use the assumption that there is a con- 
stant c > such that m^^ < cX'^j . We do not believe that this is necessary. 

2. That the result should be true is again due to the fact that the down- 
ward drift resulting from overcrowding in a site is quadratic, whereas the 
upward drift due to reproduction in the population is only linear. Moreover, 
for sufficiently crowded sites, immigration from neighboring sites is being 
compensated by intraspecific competition. 

Sketch of proof of Lemma 4.4. Again, the proof relies on a suitable 
comparison to one-dimensional diffusions. This time, the comparison is set 
up in a way such that Yi is being dominated by a one-dimensional diffusion. 
Indeed, notice that, for t > n, as long as Yi{t) > m[j/{a'\[j) for all j ^ i, an 
informal calculation shows that 

(38) dY,{t)<a'lM'-^ rn!ij-X',iYiit)\Yiit)dt + ^jY^)dBi{t). 

Hence, the immigration of mass from site j is compensated for by the down- 
ward drift due to crowding at site j and we may compare the evolution of 
Yi (again applying the Corollary A. 2 to the Ikeda-Watanabe Comparison 
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theorem) to that of the solution Zj of the logistic Feller diffusion (17), this 
time with driving Brownian motion given by {Wi{t)}t>o := {B'^{t)}t>o- It 
is then possible, in a way similar to the proof of Lemma 4.3, to construct 
the suitably measurable events E' {i,n)[v'] in terms of Zi with the required 
properties. Here, the quadratic downward drift in the logistic Feller diffusion 
does the trick, as standard speed-measure and scale-function computations 
for Zi show. Full details of this proof can again be found in Blath, Etheridge 
and Meredith [1]. □ 

Remark. Note that for a' > v'q in Lemma 4.4 (and, of course, also for 
the stronger condition a' /2 > Vq), at time n G N, recalling = 2, 

(39) Pi E'{j,n-l)[a']eT'{i,bL,n-l,2) cJ''{i,NL,n-l,2) 

implies (21). Moreover, together with the additional condition Yj{n) < a'M' 
for all + [-bL, hL^ n Z'^, both 

(40) Pi E\j,n-l)[a']eT'{i,hL,n-l,2)(ir'{i,NL,n-l,2) 
and 

(41) P E'{j,n)[a']eT'{i,hL,n,2)^T'{i,NL,n,2) 

jei+[~bL,bL]''nZ'i 

imply (21), too. Thus, all three events can be considered as instances of the 
"safe environment condition" G's^^{i,n) in the sense of Lemma 4.3. 

4.2. Comparison arguments. 

4.2.1. Comparison of C, {resp. rj) to 2N = 2{b + 2) -dependent oriented per- 
colation. We first focus on the X- (resp. (-) population and construct a 
suitable "good event" G^. Recall that by our technical assumption, L de- 
notes the maximum range of both migration matrices {mij}, {m'^} and b is 
the smallest positive integer such that the range of {jij} and {jij} is less 
than (5 — 1) • L. Also, recall from (15) the definition 

M 

if Xi{2n) > — and Yj{2n) < a'M' 

VjGz + [-feL,6L]'^nZ'^, 
otherwise 

and from the comparison assumptions (10) the notation C(n) G H if Ci{n) = 
1 for all i € [-L/2, L/2f CiZ'^ . Finally, recall from (11) the notion of trans- 
lation of H by mL, for some m G Z, denoted hy m ■ H. 



(42) 
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Lemma 4.5 (Comparison). Let G (0,1). Suppose we are given fixed 
parameters a,a' ,{mij},{m[j},{Xij},{X'^j}. Then, under the assumptions on 
Model I for Theorem 1.2, there exist finite constants Mq > and a' > 
such that if M > Mq, K := 2aMc + 1 and M' > J2j "^'ij/^' j then there is a 
finite 7 = 7(0' M') > such that ifJ2jlij < 7; ^.''^d for all {'y'ij}, the process 
{C(^)}n>o fulfils the modified comparison assumptions (14) for NL and 0. 
In particular, for each configuration Q ^ H , there exists a ^^good event" 

Gce.^*(iVL,[0,2]), 

where N = b + 2, with P[GJ >l-9, such that if C(0) = C, then on G^, 

C{i)e{+i)Hn{-i)H. 

Consequently, using Theorem 3.5, the process Af„ := {{m,n) G C:({n) € m ■ 
H} dominates a 2N -dependent oriented percolation process \yVn\n>o on C 
with density at least 1 — 6 and initial condition Wq = Xq . 



Remarks. 1. A similar result is true for the Y- (resp. rj-) popula- 
tion which, given ^ > 0, produces a similar threshold Mq and parameters 
M',a,^', which allow a comparison to a 2N = 2(6 + 2)-dependent oriented 
percolation process of density at least 1 — 9 via a similar "good event" Grj ■ 

2. The available degree of freedom in the choice of {'Jij} in this result is 
crucial for the simultaneous comparison of {Cin)}n>o and {r]{n)}n>o which 
we will need to consider later. It is due to the fact that our results for 
the "control of the environment" in Lemma 4.4 are entirely independent of 
these {7j-j}, since competition by X only facilitates the "good environment 
condition" determined in terms of Y. 



Proof of Lemma 4.5. Fix 9>0 and let 

e 



^^^^ ^ 2(4(6 + 2)L)'^- 

We begin with the specification of consistent parameter values for our model 
that will lead to the required comparison. First, note that all of the con- 
stants a,a' ,{mij},{m'^j},{Xij},{X'^j},b,N,L will remain fixed throughout 
what follows. The only values we need to adjust suitably in order to pro- 
duce the proof are M,M' ,a', {'Jij}- The proof is entirely independent of the 
choice of {7^^} (provided all parameter values remain compatible with the 
assumptions of Theorem 1.2). 

First, we choose sufficiently large a' so that for any M' > ^jm'^/a' , 

(44) a' > 



2a'M'c + 1 
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and, moreover, such that a! 12 > v'q in Lemma 4.4 with the above e. Then, 
for each i and n, we have the bound 

(45) F[E'{i,n)[a'/2]]>l-e. 

Note that this bound does not depend on {jij} (and obviously not on {'Jij})- 
From now on, a' remains fixed. Define, for any M' > m'^^/a' , the constant 
7 = 7(a', M') := (2a'M')"^ so that, for each i, we have Y.j&Z'' nfia'M' < 1. 

Finally, we can find Mq > such that for all M > Mq, the bounds of 
Lemma 4.3 for the "infection" and "nonrecovery probabilities" hold with 
our choice of e. 

We now check that with these parameter values for Mo,a',7 and for all 
M > Mo, assuming C(0) G H, there is a "good event" G T*{NL, [0,2]), 
which implies 

(46) ^(i)e(+i)//n(-i)i? 

and has probability at least 1 — 9. Recall that (^(0) € H means: 

• Xi(0) > M/K, where K = 2aMc + 1, for all i G [-L/2, L/2]'^ n Z^, 

• Yj{0)<a'M', for all j e i + [-bL,bL]'^ nZ'^,i £ [-L/2, 1/2]^^ nZ"^. 

To construct G(^, recall that one time step for C corresponds to two time 
units for X and Y. We split the corresponding time interval [0, 2] into two 
parts, [0, 1] and [1, 2]. By Lemma 4.4, applied with v' = a' > v'q, we see that 

f]E'{i,0)[a']:i e[-{b + l/2)L,{b + 1/2)1]'^ nZ" 

i 

>l-i2{b + l/2)Lfe, 

and recall that this event, denoted by ^'(0, (6+ l/2)L,0,2)[a'] G J^'(0, {b + 
1/2)L,0, 2) for short, implies, since C(0) G H, by Lemma 4.4, that 

sup {Y^(s)■.ie[-{b + l/2)L,{b + l/2)L]'^nZ'^}<2aM'. 

0<s<2 

Next, Lemma 4.3 tells us that [recall Gnonrcc (i,0) G J^(i,0,l) for n = 0]. 



nG'nonrcc(i,0) : i G [-L/2, L/2]<^ H Z« 



> 1 - L^e 



and this event, denoted by Gnonrcc(0, 1//2, 0, 1) G J^{0, L/2,0, 1), implies 



^^(1) > 



M 



for alHG [-L/2,L/2]'^n: 



2aMc+l 

Now, applying Lemma 4.4 once again, this time with v' = a' /2, we see that 



f]E'{i,0)[a'/2]:ie[-NL,NL]'^nZ'^ 



> 1 - e{2NLY 
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and this event, denoted by S'(0,iVL,0,2)[aV2] e :r'(0, A^L, 0, 2), ensures 
that 

sup {Yi{s):ie[-NL,NLfr\Z'^}<a'M'. 

sG[l,2] 

Combining all of these events, we have guaranteed that with probability at 
least 1 - {2NLfe - L'^e - (2(6 + l/2)LYe, we have that 

(47) X,(l) > ^ = ^^^^^ for all i G [-L/2, 1/2^ n Z"^ 
and 

(48) sup {Yj (t) : i e [-iVL, NL\^ n Z^} < a M' , 

l<t<2 

where N = b + 2. Thus, throughout the time interval [1,2], a "safe environ- 
ment condition" Gsgj.(i,l) G J-'{i,NL, 1,1) holds at time n=l for all i G 
[-(3/2)L, (3/2)L]'^ n and hence the local Y-population is not "too big," 
that is, for all j G [-(3/2)L, (3/2)L]'^nZ'', we have max^^ji 2] SzeZ'' 7i«^/(*) < 
1. Thus, the ("-process can safely invade the neighboring boxes, that is, con- 
ditional on (47) and given the above instance of the "safe environment con- 
dition," for each site 

(49) iG{{Lei + [-L/2,L/2]'^} U {-Lei + [-L/2, L/2]'^}} n Z^ 

where ei denotes the first unit vector in Z"^, the "infection event" Ginfec(^5 1) 
at i has probability greater than 1 — e, by Lemma 4.3 (noting that by our 
choice of L, and the fact that rriij is a function of \\i — j\\ alone, each such 
site i has at least one occupied neighbor j G [— L/2,L/2]). Hence, after all 
of these prerequisites, the probability that simultaneously for all such sites 
i taken from the set in (49), at time 1, the event 

Ginfec(«, J, 1) e ^{i. NL, 1, 1) holds for some j G [-L/2, L/2]'^ n Z'', 

implying, under the above conditions that Xi{2) > M/K = M/{2aMc + 
1), is at least 1 — 2L'^e. We denote this event by Ginfec(-^ T L/2,1,1) G 
Jf^(0,iVL,0,2). Thus, we may define the J^*(7VL, [0, 2])-measurable "good 
event" 

:= E'{0, (b + 1/2)L, 0, 2) [a'] D Gnonrcc(i, L/2, 0, 1) 

(50) 

n E'{i, NL, 0, 2) [a'/2] n Gi,,Uh L T L/2, 1,1), 
which implies, given Q = C(0) G H, that C(l) G {+1)H n { — 1)H and observe 
¥[G(^] > 1 - (2(6 + l/2)L)'^e - L'^e - {2NLfe - 2L'^e 
>l-4:{2NL)'^e = l-e, 
which completes the comparison. □ 
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4.2.2. Simultaneous comparison and proof of Theorem 1.2 and Corollary 
1.3. Assume that 9 < . We may then choose a, a' such that for 

all M > J2j rriij/oL and M' > J2j f^'ij/o-' ■, (45) holds for both populations X 
and Y with e = \ point is that this can be done simultaneously, 

since the bounds for the control of the environment do not depend on the 
behavior of the competitor. 

We may then pick M,M' and simultaneously 7,7' such that Lemma 4.3 
holds with e = jj2WLy both the X- and the F-population. Condition 
(45) is unaffected by this, since the bounds on the environment do not de- 
pend on {"fij} and {jij} and hold for all M > J2j ^ij/ct and M' > J2j ^'ij/oi' ■ 
Assuming, then, that ^(0) and r/(0) S H, observing that condition (44) on 
a' ensures that n'l := 2a'M'c+i ^ cl'M' =: 4 (with a similar inequality for 
Ki,K2) leads to the initial condition 

(x(o), y(o)) G h{k^,k2; 4,4; + 1/2)^) 

specified in Theorem 1.2, with 



M' 



K,2 



' - a'M'. 



' 2a'M'c + 1 

Hence, we can simultaneously construct the corresponding good events 
and Gn and infer from Theorem 3.2 that both the X- and Y-population 
survive, each with probability greater than which yields persistence of 
{X, y} with positive probability. Moreover, if we make the stronger assump- 
tion that the initial configurations of the X- and y-populations satisfy 

{X{0),Y{0))eH{KuK2; 4,4; oo), 

hence assuming 0(0) = 1, Vii^) — 1, for all i € Z"^, then, according to Theo- 
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rem 3.6, liminf„_^oo IP'[C2n(0) = 1] > The same result holds for rj. Thus, 



liminfP[C2n(0) = l,7?2n(0) = 1] > ^. 

By the definition of {Cn, ??n}n>o and our bounds from the last section applied 
in a similar fashion, this implies that there is a uniform positive lower bound 
on F[Xo{t) > M/{2K),Yo{t) > M/{2K) Mt G [0,4] given C2n(0) = l,r/2n(0) = 
1] and the proof is completed. 

4.2.3. Proof of Theorem 1.4. The proof of Theorem 1.4 again follows 
by comparison to suitable one-dimensional diffusions. This time, there is no 
need to control a potentially harmful "environment," making things much 
easier. For details, we again refer to Blath, Etheridge and Meredith [1]. 
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APPENDIX 



The following classical comparison theorem and, in particular, its subse- 
quent corollary, whose proof can be found in Blath, Etheridge and Meredith 
[1], are tailored for our purposes in Section 4. Note that the corollary allows 
a comparison, even if (51) only holds for intervals. 

Theorem A.l (Ikeda and Watanabe [10]). Let {Q,J',{J='t},F) be a fil- 
tered probability space and let xi{t,ij),X2{t,uj) be two real {Tt}-adapted pro- 
cesses. Let B{t,uj) be a one- dimensional {J^t}-Brownian motion such that 
B{0) = a.s. and let (5i(t,uj), (32{t,i^) be two real {Tt}-adapted previsible 
drifts. Assume that with probability one, 



and that pathwise uniqueness of solutions holds for at least one of the equa- 
tions. Moreover, assume that with probability one, 



xi(0)<X2(0), Pi{t)<bi{t,xi) and /?2 W > 62(t, ^2) Vt > 



for two real continuous functions bi{t,x),b2{t,x) on [0,oo) x M such that 



for all t>0 and x € M. Then, xi{t) < X2{t) for every t > 0. 

Corollary A. 2. In the framework of Theorem A.l, assume that xi 
and X2 are positive and nonexploding. Let 6 > 0. 

(a) Suppose that condition (51) on 61,62 is required only for all x € [5, 00). 
Assume xi(0) < ^2(0). Define t^^ := inf{t > 0:x2{t) < 6}. Then, with 
probability one, xi{t A t^^) < X2{t A rj^) for all t>0. 

(b) Suppose that condition (51) on 61,62 is required only for all x S [0,^]. 
Assume xi(0) < X2{0). Define rj^ := inf{t > 0:xi{t) > 6}. Then, with 
probability one, xi{t A t^^) <X2{t A t^^) for all t>0. 
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